DESY 04-172 ISSN 0418-9833 

MZ-TH/04-07 

|hep-ph/0410289l 

October 2004 



Inclusive Production in pp Collisions with Massive Charm 

Quarks 

B.A. Kniehl, G. Kramer, I. Schienbein 

//. Institut fur Theoretische Physik, Universitat Hamburg, 
Lumper Chaussee 149, 22761 Hamburg, Germany 

H. Spiesberger 
Institut fiir Physik, Johannes-Gutenberg-Universitat, 
Staudmger Weg 7, 55099 Mainz, Germany 

Abstract 

We calculate the next-to-leading order cross section for the inclusive production of D*^ mesons 
in pp collisions as a function of the transverse momentum and the rapidity in two approaches using 
massive or massless charm quarks. For the inclusive cross section, we derive the massless limit 
from the massive theory. We find that this limit differs from the genuine massless version with MS 
factorization by finite corrections. By adjusting subtraction terms, we establish a massive theory 
with MS subtraction which approaches the massless theory with increasing transverse momentum. 
With these results and including the contributions due to the charm and anti-charm content of 
the proton and anti-proton, we calculate the inclusive D*^ cross section in pp collisions using 
realistic evolved non-perturbative fragmentation functions and compare with recent data from the 
CDF Collaboration at the Fermilab Tevatron at center-of-mass energy yfS = 1.96 TeV. We find 
reasonable, though not perfect, agreement with the measured cross sections. 
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I. INTRODUCTION 



Recently, there has been quite some interest in the study of charm production in proton- 
anti-proton collisions at high energies, both experimentally and theoretically. The CDF 
Collaboration at the Fermilab Tevatron presented results for prompt charm meson produc- 
tion cross sections at center-of-mass energy VS = 1.96 TeV. The differential cross section 
da/dpT was measured as a function of transverse momentum (jpx) i n the central rapidity 
(y) region \y\ < 1 for inclusive production of D°, D + , D* + and D£ mesons and their charge 
conjugates For defmiteness, we shall concentrate here on D*^ mesons. However, our 
results readily carry over to any other heavy-flavored hadrons. 

On the theoretical side, various approaches for next-to-leading order (NLO) calculations 
in perturbative QCD have been applied for comparison with experimental data. In the 
so-called massless scheme , also known as zero-mass variable-flavor-number (ZM-VFN) 
scheme, which is the conventional parton model approach implemented in the modified 
minimal-subtraction (MS) scheme, the zero-mass parton approximation is applied also to 
the charm quark, although its mass m is certainly much larger than the asymptotic scale 
parameter Aqcd- in this approach, the charm quark is also an incoming parton originating 
from the proton or anti-proton, leading to additional contributions, besides those from the 
gluon g and the u, d and s quarks. The charm quark fragments into the D*^ meson similarly 
as the gluon and the light quarks with a fragmentation function (FF) known from other 
processes. The well-known factorization theorem then provides a straightforward procedure 
for order-by-order perturbative calculations. Although this approach can be used as soon 
as the factorization scales of the initial and final states are above the starting scale of the 
parton distribution functions (PDFs) of the (anti-)proton and of the FFs of the D* ± meson, 
the predictions are reliable only in the region of large transverse momenta px 3> ra, where 
terms of the order of m 2 /p T can safely be neglected. 

Another calculational scheme for heavy-flavor production which could be applied to the 
process p + p — > + X 0, 0, B] is the so-called massive scheme, also called fixed 
flavor-number (FFN) scheme, in which the number of active quark flavors in the initial state 
is limited to rif = 3 and the charm quark appears only in the final state. In this case, the 
charm quark is always treated as a heavy particle and never as a parton. The actual mass 
parameter m is explicitly taken into account along with the variable px as if they were of the 
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same order, irrespective of their actual relative magnitudes. In this scheme, the charm mass 
acts as a cutoff for the initial- and final-state collinear singularities and sets the scale for the 
perturbative calculations. However, in NLO, terms proportional to ck s (//r) lnfjo^/m 2 ), where 

is the renormalization scale, arise from collinear emissions of a gluon by the charm quark 
at large transverse momenta or from almost collinear branchings of gluons into cc pairs. 
These terms are of order unity for large values of pt and, with the choice fiR ~ pt, they 
spoil the convergence of the perturbation series. The FFN approach with rif = 3 should 
thus be limited to a rather narrow range of pt values, reaching up to a few times m. 

There are also interpolating schemes, which smoothly interpolate between the FFN 
scheme at low values of pt and the ZM-VFN scheme at large values of pt, with some 
reedom concerning the detailed implementation. The Aivazis-Collins-Olness-Tung (ACOT) 
^1 scheme, also known as general-mass variable-flavor-number (GM-VFN) schemers one of 
them. This scheme was applied to the hadroproduction of heavy flavors in Ref. [9( taking 
into account the FFN part at NLO and the leading logarithms of the ZM-VFN part. 

Another interpolating scheme which has been applied to inclusive D*^ production in the 
Tevatron region is the so-called fixed-order next-to-leading-logarithmic (FONLL) scheme, 
in which the traditional cross section in the FFN scheme and a suitably modified cross 
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section in the ZM-VFN scheme with perturbative FFs are linearly combined 
combination is done in such a way, that the ZM-VFN term is weighted with an ad-hoc 
coefficient function of the form Pt/{Pt + 25m 2 ) to enforce its suppression in the low-pr 
range. In both finite-charm-mass approaches, the FFN and the FONLL, the theoretically 
calculated FFN cross sections are convoluted with a non-perturbative FF extracted from 
e + e~ data. This assumes universality of the FF which is not supported by a factorization 
theorem as in the ZM-VFN approach. 

As has been explained at many places in the literature, mainly in the context of charm 
production in deep- inelastic ep scattering (for a recent review, see Ref. ^^), the correct 
approach for pt 3> m is to absorb the potentially large logarithms into the charm PDF of 
the (anti-)proton and the FF of the c — > D* + transition. Then, large logarithms of the type 
ln(fi^ I m 2 ) , defined with the factorization scale /i^, determine the evolution to higher scales 
and can be resummed by virtue of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
[rJ evolution equations. The unsubtracted terms of the form ln^/yU^) are of order unity 
for the appropriate choice of /i^ of order px- After factorizing the lnm 2 terms, the hard cross 
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section is infrared safe, and rif = 4 is taken in the formula for a s and the DGLAP evolution 
equations. The remaining dependence on m, i.e. the terms proportional to m 2 /p T , can be 
kept in the hard cross section to achieve better accuracy in the intermediate region px > m. 
The factorization of mass- divergent terms can be extended consistently to higher orders in 
a s , as has been shown by Collins in the context of heavy-flavor production in high-Q 2 ep 
collisions jl^ . 

It is well known that the subtraction of just the collinearly (mass) singular terms does 
not define a unique factorization prescription. Also finite terms must be specified. In 
the conventional ZM-VFN calculation, the mass m is put to zero from the beginning and 
the collinearly divergent terms are defined with the help of dimensional regularization. This 
fixes the finite terms in a specific way, and their form is inherent to the chosen regularization 
procedure. If one starts with m ^ and performs the limit m — > afterwards, the finite 
terms can be different. These terms have to be removed by subtraction together with the 
lnm 2 terms in such a way that, in the limit px — > oo, the known ZM-VFN expressions 
are recovered. This requirement is actually unavoidable, since almost all existing PDFs and 
FFs, including those for heavy flavors, are defined in this particular scheme (or sometimes in 
the deep-inelastic-scattering (DIS) scheme, which can be derived from the MS scheme). It is 
clear that a subtraction scheme defined in this way is a correct extension of the conventional 
ZM-VFN scheme to include charm-quark mass effects in a consistent way. In the following, 
we shall refer to it as the GM-VFN scheme, since it is conceptionally similar to the framework 
of Ref . P| . For a fully consistent analysis of heavy-flavor production in pp collisions, it will 
eventually be necessary to use dedicated PDFs and FFs with heavy-quark mass effects 
included, determined by global fits utilizing massive hard-scattering cross sections. Needless 
to say that it is, therefore, important to work out massive hard-scattering coefficients in 
one particular scheme for all relevant processes. Actually, just recently PDFs of the proton 
with heavy-quark mass effects included have been constructed by members of the CTEQ 
Collaboration in a scheme very similar to ours as outlined above. If these were used 
in a calculation of charm production in pp collisions, the treatment of the corresponding 
hard-scattering cross sections would have to be adjusted to these PDFs. However, we think 
that this would be premature as long as similar constructions of the FF for c — > D* + do not 
exist. 

In a recent work, two of us applied the GM-VFN scheme to the calculation of the cross 
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sections for 7 + 7 — > + X 



16 



17| and jp -»• D*± + X H In Ref. H, we considered 
only the direct and the single-resolved cross sections with In the calculation of the 

full cross section for 7 + 7 — > D*^ + X, needed for the comparison with experimental data, 
i.e. in the sum of the direct, single-resolved and double-resolved parts, the double-resolved 
contribution was still treated in the ZM-VFN scheme with rif = 4. It is the purpose of 
this work to apply the GM-VFN approach to the pp cross section. The results of this 
calculation can then also be applied to the cross sections of double-resolved 77 and resolved 
7P collisions. These cross sections play an important role due to the partonic subprocesses 
g + g —* c+c and q + q — > c+c with charm quarks in the final state and due to the subprocess 
g + <?(<?) — > c + c + q(q), where q is one of the light (massless) quarks u, d and s. These 
contributions and their NLO corrections should be computed with massive charm quarks. 
Although FFs for various charm mesons have been constructed from e + e~ data we shall 
restrict ourselves to inclusive -D* ± production and study the mass-dependent corrections for 
this special final state only. Results for the inclusive production of other charm mesons will 
be presented in a future publication. 

Starting with g + g — > c + c, the NLO corrections for this subprocess can be split into an 
Abelian and two non-Abelian parts. The Abelian part is, up to an overall constant factor, 
identical to the NLO corrections to 7 + 7 — > c + c. For this part, the terms in the massive 
theory surviving in the limit m — > 0, which are not present in the ZM-VFN approach, have 
been identified in our earlier work [lfij . Therefore, only the two non-Abelian parts of the 
NLO corrections to the gluon-gluon fusion cross section have to be investigated in addition 
to the cross sections for q + q — > c + c + g and g + q(q) — > c + c + q(q). 

The NLO corrections with non-zero quark mass m were calculated by several groups 
P Q El 2 In neither of these references, complete formulas for the NLO corrections were 
published. Fortunately, Bojak supplied us with the computer code which was used in Ref. 0. 
From this code, we were able to read off the complete NLO squared matrix elements needed 
for the computation of the mass- dependent cross section. The authors of Ref. compared 
their results with those of Refs. 0, f| and found complete agreement. Therefore, we use these 
expressions to derive the limit m — > and establish the subtraction terms by comparing 
to the MS-factorized cross section derived in Ref. j^. The latter is available to us in the 
form of a FORTRAN program Since, in the work of Ref. Bl, the FFN cross section 



was derived with a method different from the one used in Ref. 
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space slicing method for separating the infrared-divergent part from the hard part of the 
cross section, we also derive the massless limit of the Abelian part for consistency. With 
this knowledge, we can compute the finite-mass corrections for the full NLO cross section 
with MS factorization. 

The outline of our work is as follows. In Sec. |HJ we describe the formulae which we use 
to calculate the cross section for g + g — > c(c) + X, q + q — > c(c) + X, g + q — > c(c) + X and 
g + q — > c(c) + X with non-zero charm-quark mass. For these cross sections, we perform the 
limit m — > and compare the results with the ZM-VFN theory of Ref. j^J. The results are 
collected in Sec. IIHI and three appendices. In Sec. we also present numerical results to 
test the validity of the subtraction terms and show how the various terms in the NLO cross 
section approach their corresponding massless limits for large values of px- After adding the 
contributions with (anti-) charm quarks in the initial state, which are present in the ZM-VFN 
scheme with nt = 4, as well as the contributions due to the fragmentation of gluons and 

n 

light (anti-) quarks, we compare our results to recent experimental data from CDF |1| in 
Sec. IIVI A summary and conclusions are given in Sec. 

II. LO AND NLO DIFFERENTIAL CROSS SECTIONS 

The differential inclusive cross section for the process p + p — > D*^ + X has many con- 
tributions. In this section, we consider those contributions where the charm quark appears 
only in the final state. We study the charm-quark mass dependence to obtain the massless 
limit, which is then compared with the ZM-VFN theory, and to establish the influence of 
the m 2 /p T terms in the GM-VFN theory defined in the same MS factorization scheme as 
the ZM-VFN theory. 

There are only two leading-order (LO) partonic subprocesses, g + g — > c + c and q + q —* 
c + c. The NLO corrections to these two channels comprise the virtual corrections and 
gluonic bremsstrahlung contributions, g + g — > c + c + g and q + q—^c + c + g. In addition, 
the subprocesses g + q — > c + c + q and g + q^c + c + q appear for the first time at NLO. 
In the following subsections, we present the LO cross sections in order to fix the notation. 
Then, we explain how we calculate the NLO corrections to g + g c + c and q + q —>■ c + c 
and the cross sections for g + q(q) — >• c + c + q(q). 
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A. LO Cross Section 



We start with the subprocess 

g(h) + 9(h) -> c(pi) + c(p 2 ) + [^(p 3 )] , (l) 

where ki, k 2 and (i = 1,2,3) denote the four-momenta of the incoming gluons, the 
outgoing charm and anti-charm quarks and a possible gluon in the final state (in square 
brackets). We have the following invariants 

s = (h + k 2 f , 

ti — t — m 2 = (ki — pi) 2 — m 2 , 

Ui = u — m 2 = (k 2 — Pi) 2 — m 2 , 

S2 = (h + k 2 - pi) 2 - m 2 = s + ti + ui . (2) 

Here ti and ui are determined by the four-momentum of the observed charm quark. As 
usual, we define the dimensionless variables v and w by 

v = l + —, w = — , (3) 

S S + ti 

so that ti = —s(l — v), Ui = —svw and s 2 = sv(l — w). For p 3 = 0, i.e. at LO, we have 
s 2 = and w — 1. 

The LO cross section for g + g ^ c + c is 



d Of o . \ r / -i X I ^ ^ hui 

where 



U Ui Am s ( m s 
— + — + 1 

Ui ti tiUi \ tiUi 



(4) 
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and all color factors have been expressed in terms of the Casimir operators Cf = (N 2 — 
1)/(2N) and Ca = N, where denotes the number of colors. 
The LO cross section for q + q — > c + c reads 

=«,(.),(! -„) C ^A_L + _j, (6) 

where 
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B. NLO Cross Section 

Following the notation of Ref. Q, the color decomposition of the NLO squared matrix 
element for g + g — > c + c + g can be written as 



M 99 \ 2 = g«E? 



N 2 - 1 



Cf-Dqed + -C a D Q q + -Ca {Ca — 2Cp) -Dkq 



(8) 



where g 2 = Ana s is the strong coupling and E e = 1/(1 — e) originates from averaging over 
the gluon spins in n = 4 — 2e space-time dimensions. The squared matrix element for the 
virtual corrections to g + g — > c + c can be written in a similar fashion, there is, however, 
an additional contribution from quark loops which comes with a color factor The 
Abelian contribution -Dqed is identical to the QED part of 75 — > ccg in Ref. 22]. In 
addition, we have two non-Abelian oarts D OQ and D KQ . For isolating the divergences in 
the soft limit, Bojak and Stratmann [2J used the same method as in Ref. p|. They slice the 
2^3 contributions into a soft- gluon and a hard-gluon part by introducing a small auxiliary 
quantity A. In the limit A — ► 0, the kinematics of the soft-gluon cross section is that 
of the 2 — > 2 process, so that the phase-space integrations can be performed analytically. 
After combination with the virtual cross section, the infrared 1/e and the combined infrared- 
collinear 1/e 2 singularities, all proportional to the n-dimensional LO cross section, cancel. In 
this way, the soft plus virtual cross section becomes finite, except for the remaining collinear 
1/e singularities, which cancel against the subtraction terms in the collinear factorization 
procedure of the gluon PDF of the (anti-)proton. The integration over the phase space of 
the two unobserved partons in the hard part is done analytically as far as possible, using 
the methods of Refs. 0, f|. 

The same steps are taken to calculate the NLO corrections to the subprocess q+q — > c+c. 
The squared matrix element for the real corrections is color-decomposed in the following way 

\M qq \ 2 = g 6 ^ (c 2 F N QED + \c F C A N ^j , (9) 

where iVqEp is again the Abelian and iVoK the non-Abelian part, which are also obtained 
from Ref. 7(. A similar decomposition is used for the virtual corrections, which receive 
an additional contribution from one diagram with a quark loop in the gluon propagator, 
proportional to the color factor Cf/2. 
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The subprocess 

g(h) + q{h) -> c(pi) + c(pa) + g(p 3 ) , (10) 

which occurs only in NLO, has two pieces in color space. The squared matrix element is 
split up according to 

\M gg \ 2 = g 6 E t ^ (c F J QED + \C A J ^\ . (11) 

The Abelian part Jqed is, up to a factor, equal to the corresponding squared matrix element 
for the subprocess 7 + <? — > c + c + q. The squared matrix element for the crossed subprocess 
g + q — > c + c + q, i.e. Eq. (fTU|) with the quark replaced by an anti-quark, has the same 
structure as in Eq. (fTTj). but with slightly different coefficients. 



III. ZERO-MASS LIMIT OF THE MASSIVE CROSS SECTIONS 



In this section, we collect our results for the cross sections in the limit m — > 0. We 
consider the following contributions: (i) the NLO corrections to g + g — > c + c, (ii) the NLO 
corrections to g + g— >c + c and (iii) the process g + q-^ c + c + q and the corresponding 
channel with q —>■ q, where q denotes any of the light (massless) quarks u, d and s. The 
result for the limit m — > will in general be different from the cross section obtained in the 
ZM-VFN approach, where the mass of the charm quark is neglected from the beginning. In 
the genuine ZM-VFN calculation, originally performed by Aversa et al. 20j, the collinear 
singularities connected with the charm quark appear as 1/e poles in dimensional regular- 
ization. In the FFN theory, they appear as terms proportional to ln(m 2 /s), instead. So, in 
this theory, the collinear singularities are regularized with a finite charm mass. Due to this 
different procedure for regularizing the collinearly divergent contributions, different finite 
terms appear. The occurrence of different finite terms in these two regularization schemes 
is due to the fact that the two limits, m — > and e — > 0, are not interchangeable. If one 
wants to implement the factorization of these collinearly singular terms in the MS scheme, 
as is done in the ZM-VFN scheme from the start, the different finite terms, which turn up 
in the limit m — > 0, must be subtracted. Such finite terms have already been found for the 

n 

case of the NLO corrections to7 + 7^c + c |16j and 7 + g — > c + c and for the cross section 
oi^j + q-^c + c + q |l7|. 



9 



As in Refs. [l^L fl7 1 . we decompose the NLO contributions to the cross section in the limit 
m —>■ as follows: 

lim — — - — = I ci + c\ In — 5(1 — w) 

m-*o dvdw \ s J 

+ c 2 + c 2 In — + c 3 ' 



+ c 5 In v + Cq ln(l — vw) + C7 ln(l — v + vw) + c 8 ln(l 



2 

m 



+ eg In w + cio ln(l — w) + en + en In 

ln(l-i7 + t;w) lnw ln (i^w) no s 

+ C12 I-C13- hcu— • (12) 

1 — w 1 — w 1 — w 

The coefficients q are mass-independent; the dependence on the heavy-quark mass enters 
only through the logarithms ln(m 2 /s) and the choice [Ir = fip = m. 

In the following subsections, we shall present the complete expressions for the coefficients 
q. The massless limit of the cross sections from Ref. |7| will be compared with the results 
of Ref. Q obtained in the ZM-VFN theory. 



A. Massless Limit of NLO Corrections to g -\- g — ►c + c 

Before we write down the coefficients q, we give the LO cross section for g + g 
with m — 0. It has the simple form 

d 2 af 9 



'LO 



m->o dvdw 



c(s)5(l - w)t(v)[C f - C A v(l - v)} 



with 



t(v) 



v 1 — V 

+ 



1-v 



The results for the various coefficients q are written in the form 



c + c 



(13) 



(14) 



(15) 



where o, are the results of Ref. 201 in the ZM-VFN scheme and Ac,- are the subtraction terms 
needed to convert the cross section of Ref. Q to the GM-VFN scheme. The coefficients are 
decomposed into one Abelian part, two non-Abelian parts and a quark- loop contribution 
(which only occurs for C\) in the following way: 



Ci = C{s) 



Cl cT d + \C\ cT + \C A (C A - 2CV) + 5 a ^cf 



(16) 
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with 



C(s) 



a: 



-c(s 



(17) 



2(N 2 - l)s 2tt 

The expressions for the q are lengthy and, therefore, delegated to Appendix 1X1 

The coefficients q, which agree with results obtained from Ref. [20] , refer to the version 
where every incoming gluon is averaged with the factor 1/[2(1 — e)] and where the factor- 
ization of singularities due to collinear quarks and gluons is p erformed in the customary 
MS subtraction scheme. Furthermore, deviating from Ref. 20(, in the expressions for Cj in 
Appendix El the renormalization scale fin and the initial- and final-state factorization scales 
fip and fi' F are identified with m, i.e. = /ip = fi' F = m. The non-vanishing subtraction 
terms in the massless limit of the FFN theory of Ref. |?J are found in Appendix EI They are 
Aci, Ac 2 , Ac 3 , Ac 5 , Acio an d Acn. For all three contributions c qed , c oq and cj" 1 , we have 
the relation 

Ac 5 = Ac 10 = 2Acn . (18) 



(19) 



Denoting the ZM-NLO result of Ref. |20| by d azM/dvdw, we can thus write 



n . d ctnlo 
hm — — - — 

m-+o dvdw 



d 2 (TzM 

dvdw 



fi F = flp = Tflj + 



^ 2 0"sub 

dvdw ' 



where 



d 2 cr svl Y> 
dvdw 



Aci5(l -w) + Ac 2 



w 



Ac 3 



ln(l 



1 — w 



Ac, 



lnt> + ln(l 



w 



(20) 



For the first three subtraction terms proportional to Aci, Ac 2 and AC3, one obtains simple 
expressions, if the three contributions in Eq. (fT^j) proportional to the color factors C F , C\ 
and Ca {Ca — 2C F ) are added. From the results in Appendix |XJ we obtain 



Aci = (l-lnw-ln 2 *;) x 2C(s)C f t(v)[C f -C a v(1-v)} 
Ac 2 = -(21nu + l) x 2C(s)C F r(?;)[C i r -C>(1 - u)] 
Ac 3 = -2 x 2C(s)C f t(v)[C f - C a v(1 - v)] 



(21) 
(22) 
(23) 



and 



Ac 5 = C(s)Cj, ( C F Acf d + -C A A^ 



oq 



(24) 
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In the last equation, we have used Ac 5 q = — Acg q ; the explicit expressions for these coeffi- 
cients are given in Appendices IA 21 and IA31 Note that the last factors in Eqs. (|2T| - (j23|) are 
proportional to the LO cross section. Finally, we have to subtract the quark-loop contribu- 
tion, which is absent in the ZM-VFN scheme, via j^] 

A Cl Aci + C(s)—Acf = A Cl - C(s)C A -v(l - v) . (25) 
4 9 

The FFN theory for g + g c + X in the limit m — > approaches the ZM-VFN theory 
with scales fiR = fip = = 171 if the finite terms Ad, Ac2, AC3 and Acs = Ado = 2Acn 
as given in Eqs. ([2T |l -(|2^ jl are subtracted. As already mentioned above, the necessity for 
such a subtraction is to be expected, since the regularization of collinear singularities with a 
mass parameter m does not give the same result as the one with dimensional regularization 
and m = from the start. 

In Ref. jit] ], it was shown that the finite subtraction terms for the subprocess 7+7 — > c+X 
can be obtained by a convolution of the LO cross section with a perturbative partonic FF 
(fcjx, jj) for the transition from a massless to a massive charm quark of the following form 



d c Jx, jj) = C 



a s [ 1 + x 



2 



2tt \ 1 - x I m 2 v ; J J + v ; 

Therefore, all Abelian terms proportional to in Eqs. (|21|) - (j24|) can also be generated 
in this way. Our explicit calculations show that also all terms proportional to CfCa in 
the equations above can be obtained as final-state interaction contributions with d^(x, /j,) in 
Eq. ()26|) (without the term proportional to ln(/z 2 /m 2 ), of course). 

It is important to understand that the terms containing logarithms ln(m 2 /s), i.e. the 
coefficients c\, c<i and cii, have two different origins. On the one hand, these terms are due 
to on-shell internal charm-quark lines in the Feynman diagrams which become singular for 
m —>■ 0. Also internal gluon and light-quark lines give rise to singular contributions for m — > 
0, which have to be factorized into corresponding perturbative FFs describing the transition 
from a gluon or a light quark to the heavy charm quark. In the ZM-VFN calculation, these 
contributions have been factorized as final-state singularities and are recovered by setting 
fi' F = m. Another part of these singular terms can be assigned to the initial state and is 
found in the ZM-VFN calculation if one sets fip = m. On the other hand, in both the FFN 
and ZM-VFN calculations, there are singularities due to internal gluon lines which are, in 
both approaches, factorized as initial-state singularities. Finally, there are logarithms due 
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to the renormalization of a s . In Appendix El these terms are written down for the choice 
I^r = Hf = Hf = m - 

For our application, choosing the scales equal to the heavy-quark mass is not appropriate. 
For the case of large px values, a common choice is [Ir = fip = fi' F = £my, where my = 
a/ m 2 + p T is the transverse mass of the meson and £ = 0(1). Therefore, we have to 
rescale all terms proportional to lnm 2 . In the following, we shall give the necessary terms 
for the conversion to arbitrary scales. 

First, we present the contributions related to renormalization and initial-state factoriza- 
tion of singularities due to internal on-shell gluon lines. These terms are present in the FFN 
calculation and can, therefore, be obtained from Ref. Q in the limit m — > 0. The rescaling 
is obtained by adding the following terms to the cross section in Eq. Q12j) : 

(^f) = ( k\n^ 2 +^1^ + 1^^) 5(1 -w) 

\dvdwj rcscal \ m 2 m 2 m 2 J 

+ (d 2 \n^ + l\n^)(-^—) +d u ln4 + liln4- ( 27 ) 

The non-zero coefficients read 

dx = C(s)r(v)2pi n, - 1) [C F - C A v(l - v)} , (28) 
d x = C(s)t(v) [C f - C A v(l - v)} [2C A [\n(l - v) - ln(«)] - 2^ n/_1) } , (29) 
d 2 = -C(s)t(v)4C a [C f - C A v(l - v)} , (30) 

C\ Toq , C A (C A — 2C F 



dn = C(s 
with I3^ f) = 11N/6 - n//3 and 



Joq , u i K^A - ^F) Jkq 
4 11 4 11 



~ oq A(-l + hv - 12v 2 + 2®v 3 - 12v 4 + Av 5 ) 
dn — 2 

vfv 

4(1 - Av + 8v 2 - 8v 3 + Av 4 ) 4(3 - 9v + Uv 2 - llv 3 + Av 4 



Vivw 2 vfw 
4(-l + 5v- 12v 2 + 28v 3 - 2Av 4 + 12v b )w lQv 3 (l + v)w 2 



+ 



(31) 



2 1 2 

vfv vf 
16v 4 w 3 Av x v Av(-A + 3v) Av (10 - 13v + 5v 2 ) 

~^l~ Jr lo + x~ 2 + vfx ' {d2) 

~k q -4(-l + 3v-Av 2 + Av 3 ) 4(1 - 2v + 2v 2 ) 4(1 - v + v 3 ) 

d lx — 2 I 5 ' 2 

VfV V\VW Z vfw 

4(-l + 3v - Av 2 + Av 3 )w Av x v Av 2 v Av(A - 3v + v 2 ) 
+ vfv X^+X 2 " v\~X ' (33) 
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where we have used the abbreviations X = 1 — vw and Vi = i — v. We repeat that the 
expressions above are given in the limit m — > and will be used in our ZM-VFN calculation 
with rescaling. In the FFN calculation, we shall use instead the corresponding terms as 



given in lief. p| including their full mass dependence. 

Finally, the remaining logarithms emerging in the limit m — > due to internal charm- 
quark lines becoming massless are again found by comparing the massless limit of the FFN 
cross section with the ZM-VFN cross section. These terms are associated with a rescaling 
of the subtraction terms in Eq. ([20)1. and, therefore, we write them in the following form: 



(tT-1 = fArf 1 ln4 + Arfiln4 + A ^ ln 4) 

\dvdw J Arescal V m2 m2 ™> 2 J 



6(1- w) 



+ ( Ad 2 In 4 + Al 2 In 4) ( T~ — 1 + Ad n In 4 + AJ U In 4 • (34) 

These terms must be added to Eq. (|2()jl . i.e. subtracted from Eq. p2jl . to rescale to the 
appropriate renormalization and factorization scales if one wants to use the FF for the 
transition c — ► D* + and the charm PDF of the (anti-)proton. Since, at present, we have at 
our disposal only PDFs and FFs, which are based on a ZM-VFN calculation, we shall take 
the corresponding coefficients for m = as well. Note that this does not entail any loss 
of accuracy, as has been discussed in Ref. j3] in the context of deep-inelastic scattering. 
Moreover, this fact is of great practical importance, since the known coefficients of the ZM- 
VFN scheme, e.g. those of Ref. [20(, can simply be used, whereas their massive counterparts 
are unknown and can only be obtained through a dedicated calculation. The non-zero 
coefficients are given by 

Adi = \c{s)t{v) [C f - C A v(l - v)} , (35) 

Adi = -\c{s)t{v) [C f - C>(1 - v)\ , (36) 

Ad, = C(s)t{v)C f (^2\nv + [C F - C A v(l - v)} , (37) 

Ad 2 = 2C{s)t(v)C f [C f - C A v{l - v)\ , (38) 
Ad n = C{s) 



d Adf* + S Ad\l + ° A {Ca ~ 2Cf) Ad$ 



Ad u = C(s) 



d Adjf + °± Adg + ° A {Ca ~ 2Cf) Aft 



(39) 
(40) 
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with 



A ~ d = _ L + 2-2v + v 2 _v 2 w_2_v^ 

V\ vw vi Y 

a fa on OA \ ■ 8^i(l-2^ + 2t; 2 ) 16(l-3z; + 2z; 2 ) 
Ad^ = 2u(l - 24w) H h 



8w 2 (7- 14f + 8f 2 )w 16v 3 (-l + 2z;)w 2 16v 4 w 3 



+ " o ' " 5 + 



2 1 2 

Avv\ _ - Hv + 5v 2 ) 2^(25 - 18v + 4w 2 ) 
+ ^3 y~2 ' Y ' ' ' 

A j5 = _ 2 „ + _ + 2 "' 3 - 6 ; +4 " 2 » , (43) 

a Toed 2 - 2z; + 3v 2 - 4v 3 1 + v 2 (-2 + 2v-2v 2 + 3v 3 )w 
Ad?? = 1 — 

flf V\V 

2vvi 2v v(-3 + 4v-2v 2 ) 
2(2-z; + 3t> 2 ) 4(l-2t> + 2z; 2 ) 



Ad^ = — > ^ + 



V\VW 2 



2(-3 + 10v - I3v 2 + 4v 3 ) 4(l + t> 2 W 4v x v 4v(-l + 2v) , . 

2 + 2 + ~F2 V ' ( 45 ) 

v(w v{ X 1 X 

Adg = -AdJ?, (46) 

where Y = 1 — v + mu. 

Comparing with the results for 7 + 7 — » c + X l lq , we see that the Abelian parts of 
Eqs. (J37)) and (j%B} agree with — ci and —£2 in Ref. jlf| (if the C(s) factors are set to unity). 
Furthermore, Adf[ d agrees with — en and Ad^ with — c n in Ref. [r|, as one would expect. 
The minus sign is due to the different convention for the logarithms in Eq. (42) in Ref. 
as compared to Eq. (J54*j) . We note that Ad^ d = — Ac^ d and Ad 1 ^ = — Ac^ (see Appendices 
IA II and I A 31 respectively). Finally, we emphasize again that the rescaling of the lnm 2 terms 
to arbitrary scales //r, //p and /z' F is achieved by adding Eq. (}2Tj) to and subtracting Eq. (|54*)l 
from the cross section in Eq. (|12j) . where all coefficients of Eqs. (I28j) - (j4(jj) have to be taken 
into account. Note also that Adi and Ad% cancel if /z# = /zp. 



B. Massless Limit of NLO Corrections toq + q — ► c + c 

In this section, we give the results for the massless limit of the hard-scattering cross 
sections for the NLO corrections to the subprocess q + q — > c + c, where q is any of the light 
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quarks assumed to be massless. The LO cross section of this process for to = is 

d 2 a™ 

where 

r q (v) = (1 - vf + v 2 . (48) 

The NLO cross section is again decomposed as indicated in Eq. (|T2|) with coefficients q 
which receive contributions from two color factors and from virtual-quark loops (which only 
appear for c\ and £i). Specifically, we have 

Ci = C g (s)^- (c F cf + C A cT + S a cf) , (49) 

where 

°'M = J; = (50) 

The expressions for cf, c^ a and cf may be found in Appendix [Bj According to these results, 
finite subtraction terms in the massless limit of Ref. [7( are present in c\, C2, C3, C5, C10 and 
Cn, and we have Ac 5 = Ac w = 2Acn as in the preceding subsection. The first three Aq 
terms have the following simple form: 

Aci = (l-\nv -\n 2 v) x 2C 9 (s)r 9 (i;)C| , (51) 

Ac 2 = -(l + 2\nv) x 2C7,(s)r,(v)C|. , (52) 

Ac 3 = -2 x 2C q (s)r q (v)C 2 F , (53) 

and 

^ , ^1 ( 2vv 2 2v 2 v 1 3v-6v 2 + 4:v 3 \ 

Ac 5 = 2C q (s)C 2 F [v-^ + ^ J . (54) 

The complete subtraction contribution for the qq cross section has the form of Eq. (|2T)|) with 
the Acj terms given in Eqs. (|51| ) -([54 p and Appendix El Again, the subtraction terms agree 
with the convolution of the perturbative FF of Eq. (|26j) with the LO qq cross section. We 
remark that the Aq terms only appear in the QED part. 

As in the gg channel, the results in Appendix [B] are given for the choice = fip = 
fj,' F = to. For our application, we need the cross section for arbitrary scales. To make the 
transition, we use Eq. (|2*7jl . where the non-zero contributions to the coefficients di, d~i and 
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di have the following form: 



d\ 
di 

du 



(«/-!) 



-C q (s)r q (y)C%[3 + 2 In v - 2 ln(l - v)] 



-AC q (s)r q (v)C 2 , 
C q (s)C^- 8V+12v2 



6ir 



l-2v + 2v< 



+ 



2V2V 2 W 
Vl 



Vl 

2v 3 w 2 v 
vi + X 



w 



(55) 
(56) 
(57) 

(58) 



Of course, for the calculation of the FFN cross section, the corresponding mass- dependent 
contributions are used as given in Ref. ij. These terms rescale all the lnm 2 contributions 
due to the renormalization and the factorization of internal light-quark lines. Therefore, the 
shift from renormalization due to internal quark loops is proportional to n/ — 1. 

To eliminate all In m 2 terms, we must also take into account the contribution from charm- 
quark loops and the one from internal charm-quark lines that is factorized into the FF. These 
terms are taken into account by using Eq. (|34j). and the corresponding coefficients are given 
by 



where 



'cf 



Ad\ 



Ad™ 



Ada 

Ada 
Ad 2 



Cq(s)T q (v)-C F , 

C q (s)r q (v)^(3 + 4lnv) 
C q (s)r q (v)2C 2 F , 



Ad 



li 



C q (s) 



2 11 



C f Ca 



Ad™ 



2v(2 - 3v + 2v 2 ) 2(1 - 2t> + 2v 2 ) 



Av 3 w Av 3 w 2 
+ 



+ 



4vvf 



Vi W 

2 Av 2 vt 2v(l - 6v + 4v 2 ) 



Vl 



Vl 



Y 2 



Y 



„ „ 4vv 2 4v(3 - 5f + 2v 2 
4v 2 v - 4v z w + 



2u(9- 12f + Av 2 



(59) 

(60) 
(61) 
(62) 



(63) 
(64) 



Y% y 2 Y 

In the general case, the contributions written down in Eqs. ()55j ) -([58 p have to be added and 
the contributions in Eqs. (|39 ^1 -(|52* ]l have to be subtracted. 
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C. Massless Limit of g + q — ► c + c + g and g + q — ►c + c + q 

The processes g + q — > c + c + q and g + q^c + c + q appear for the first time at 
NLO. In the massless limit, the corresponding cross sections are decomposed as explained 
in Eq. (|12|). Since there is no LO part, one has c\ = c<i = C3 = 0. The remaining coefficients 
are decomposed with respect to the color factors Cf and Ca as 

Cl = C q {s){C F cf + C A cT) . (65) 

The resulting coefficients cf and cf 1 are given in Appendix for the case that the observed 
transverse momentum is due to the charm quark. If the anti-charm quark is observed with 
a given value of pp, the coefficients are different. The differences c c ~ c = Ci(c observed) — 
Cj(c observed) are listed in Appendix [DJ Charge conjugation invariance relates the cross 
sections for quarks and anti-quarks in the initial state as da/dpp(g + q — > c b se rved + c + q) = 
da/dp T (g + q -> c obse rv e d + c + q). 

We find that, for all gq and gq channels, there are no extra finite subtraction terms, i.e. 
we have Aq = in all cases. Comparing the coefficients cf and cf' c ~ c with the results for 
7 + g— >c + c + g given in Ref. (l7|, we find agreement, except for Aci 1 , which was found 
to be non-zero for 7 + q — > c + c + g. The subtraction terms in Ref. II 71 were obtained by 
matching the FFN result of Ref. [2^ with the ZM-VFN result of Ref. bfil . In order to clarify 
this mismatch, we compare the formula for 7 + 0— ^c + c + o of Ref. 126] with the one from 
Ref. [27], which we extracted from the FORTRAN program used in that paper, to find that 



the two disagree. On the other hand, the result in Ref. [27J agrees with the Abelian part of 
the result for g + q — > c + c + q in Ref. (^(J. Relying on Ref. 27], we then conclude that we 



also have Acn = 0for7 + g^c + c + gin Ref. [17 1. 

As in the previous cases, the formulas written down in Appendix O are for the scale 
choice (jlr = hf = h'r = rn. The transition to arbitrary scales /zr, \ip and \il F is obtained, 
as above, using Eq. (|?7jl with the appropriate coefficients. In the present case, d\ and d 2 
do not contribute, only dn is non-zero. If we factor out C q (s), Cp and Ca as in Eq. 
we obtain the coefficients df and d™, which correspond to the initial-state factorization of 



18 



singularities due to internal gluons, in the following form: 

~ f 1 - Av + 9v 2 - Qv 3 + 2v 4 -l + 2v-Av 2 + 3v 3 -v 4 

d? - 



11 ~~ .,2 ^ „,2 



V\W 



1 - Av + 9v 2 - Qv 3 + 2v 4 )w v v(3 - v) 

^ o v ' \ ) 



j,.,, v 2 v 2 (2 + v 2 )w 2v 3 (l + v)w 2 2v 4 w" v ,_ 



v\ 2X 2 2v x X 

J 4 W 3 

v'f v'f v'f vf ' 2X 

In addition, we need the rescaling terms associated with initial-state singularities of in- 
ternal charm-quark lines and of final-state singularities of gluon lines splitting into cc pairs. 
These contributions vanish for the choice = n' F = m, which is used in Appendix O To 
convert the cross sections to arbitrary factorization scales //p and fi' F , we have to take into 
account the additional rescaling terms denoted as Adi and Adi, 

Adg = \^{l-2w + 2w 2 ) 1 (68) 
a irf 9/ \ l-2v + 2v 2 viv v(3-2v) ,„ n . 

aJ = _ v2 ^(2-4^> _ M-1 + 2^K 2«V ^ 



D. Numerical Test of Subtraction Terms and Study of Mass-Dependent Correc- 
tions 

The calculation of the subtraction terms, in particular those in Subsections IIII Al and 
IIII Bl was rather involved. Special care had to be exercised in order to recover all the 
terms proportional to 5(1 — to), 1/(1 — u>)+, [ln(l — w)j{\ — w)] + and the remaining terms 
in the decomposition of Eq. (|12jl. In particular, the delta- function terms were difficult to 
calculate, since they received contributions from several places, the virtual corrections and 
both the soft and hard parts of the real corrections. Some of these contributions contained 
a dependence on the slicing parameter. The cancellation of these contributions in the final 
results provided a partial check of our analytical calculations, in particular for the coefficients 
of the two plus-distributions. 

In order to check that the Aq terms presented in Subsections IIII Al and IIII Bl are correct, 
and also to see how the various contributions to d 2 a ab / dvdw (where ab stands for the channels 
99 1 QQi 9Q an d gq) written down in Subsection III Bl behave as functions of the transverse 
momentum p F of the D*^ meson, we calculate the NLO corrections in four different ways: 
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1) Using the results of Ref. Q], we calculate the cross sections for the channels gg (g+g — > 
c + X), qq (q + q — > c + X) and gq (g + q(q) — > c + X) in the FFN scheme. The 
corresponding results will be labeled by a m (massive calculation). 

2) We calculate the same cross sections in the limit m —>■ 0, using Eq. (|T2^) with the 
corresponding coefficients as given in Subsection IIHI and the appendices. Notice that 
m is kept at its physical value whenever it appears logarithmically in Eq. 1)12)1 . The 
results will be denoted gq (massless limit of the massive calculation). 

For both the massive calculation and its massless limit, we shall then consider the effect of 
the subtraction terms described by the coefficients Acf 

3) The subtracted massive calculation, denoted by <r^, is obtained by subtracting Eq. (|2U)l 
with the corresponding coefficients written down in the appendices from a m . This 
corresponds to the GM-VFN scheme. 

4) The subtracted massless cross sections are calculated by subtracting Eq. (|2TJ)l from <7o 
and will be denoted by a^. This prescription is identical to the ZM-VFN scheme. 

In all cases, we start from a calculation using /ir = fip = yJ F = m and take into account the 
rescaling to other renormalization and factorization scales in two steps: 

• Firstly, we rescale to //r = \lf = n'p = 2rriT adding the terms due to renormalization 
and initial-state factorization of singularities related to internal gluon lines, according 
to Eq. (|27|). This scale choice is to prevent the value of \J F from falling below the 
starting scale 2m for the fj,' F evolution of the FF for low values of pt- 



• Finally, we subtract the remaining rescaling terms of Eq. ()34j) related to internal 
charm-quark lines which become singular in the limit m —>■ 0. 

In the following, we always normalize the cross sections to the LO ZM-VFN cross sections 
calculated from Eqs. (jD) and ©, i.e. we consider the cross section ratios R a b given by 

da ab /dp T 

ab dai b o /dp T (m = 0)- { ' 

For the gq and gq channels, we normalize to the LO ZM-VFN gg cross section, since, for the 
gq channel, there is no LO contribution. All results are given at the hadron level, i.e. the 
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partonic cross sections are convoluted with the (anti-)proton PDFs and the FF for c — > D* + . 
We average over D* + and D*~ mesons. We use the PDF set CTEQ6M Q and the FF set 
obtained in Ref. Q by fitting OPAL data at NLO. Although the CTEQ6M PDFs were 
determined in the ZM-VFN scheme with rif = 5, for the time being, we only include three 
light quark flavors along with the gluon as incoming partons for all values of fip. Since, at this 
point, we wish to focus on effects related to the hard-scattering cross sections, the LO cross 
sections are evaluated using the same conventions concerning a s (/ZR), proton PDFs and D*^ 
FF as in NLO. We consider da/dp T at y/S = 1.96 TeV function of pr with y integrated 
over the range —1.0 < y < 1.0. As for the QCD input parameters, we take the charm- 
quark mass to be m = 1.5 GeV and evaluate aT f (hr) with n/ = 4 and asymptotic scale 
parameter Aq^, d = 328 MeV, which corresponds to ai (mz) = 0.1181. The renormalization 
scale fiR and the factorization scales /xp and fi' F are set equal, /i# = /ip = fi' F . Details for the 
calculation of da/dpx from d 2 a/dvdw have been given for direct 77 scattering in Eq. (45) 
of Ref. 16]. In this equation, the photon distribution functions must be replaced by the 
(anti-)proton PDFs of the gluon and the light (anti-) quarks. 

The results for the gg channel are shown in Fig. ^ In Fig. [T^, the renormalization 
and factorization scales are /Jr — fi F — ij! f — m, but the scales in a s , the PDFs and the 
FF are fixed at 2.1m, to stay above the starting scale of the FF. The full line in Fig. ^1 
corresponds to the massless limit of the FFN calculation as derived in Sec. IIII Al (crn). and 
the dashed curve is the result of the massive calculation (<r m ). We see that the massive cross 
section approaches the massless one very slowly at large values of px- At pr = 20 GeV, the 
difference between the massive and massless cross sections is still of the order of 6%. The 
ratio R gg for the massive cross section is always larger than its massless limit in the p? range 
between 5 and 100 GeV. From this comparison, we conclude that, in the gg channel, the 
terms proportional to m 2 jp\ are particularly large and lead to an increase of the massive 
cross section as compared to the massless approximation. Similar observations were made 
in Ref. Q|, where the massive and massless calculations were compared as functions of the 
mass m for fixed value of px- For our application, we are interested in the massive and 
massless cross sections, where the finite Acj terms derived in Sec. IIII Al are subtracted. This 
leads to the dashed-dotted (cr^) and the dotted (a^) curves in Fig.[T^L. We have checked that 
our result for the massless calculation after subtraction, i.e. the dotted curve, is in perfect 
agreement with the results in the ZM-VFN scheme obtained using the FORTRAN program of 
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FIG. 1: gg contribution to pp — ► Z?* ± + X including Abelian and non-Abelian parts, normalized 
to the LO cross section with m = 0. (a): Renormalization and factorization scales are /j,r = 
[if = fi'p = m (but fixed at 2.1m in a s , PDFs and FF). (b): Same as in part (a), but for 
= li F = fJ-'p = 2i7if including rescaling due to renormalization and initial-state factorization of 
singularities originating from internal gluons. (c): Same as in part (a), but including full rescaling 
to = hf = n' F = 2m.T- 

Refs. . This comparison demonstrates that the finite Aq terms, which, if subtracted, 

should produce the ZM-VFN cross section, are correct. Their subtraction from the FFN 
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result will give the GM-VFN result, which approaches the ZM-VFN result at large values 
of pt- We see that the contribution of these finite terms is by no means negligible. 

The ratios plotted in Fig. show that, in the low- to medium-pr range, the NLO cross 
section is up to a factor of about five larger than the LO one (note that the numerator of 
R gg is the sum of the LO result and the NLO corrections). This is not surprising, since we 
have chosen very low values for the renormalization and factorization scales. Since we are 
interested in the region where pr ^ Tn, a better choice of scales is fi ~ px- As usual, we 
choose hr = hf = Hp = ^ m T, which can be used for both small and large values of p?. To 
obtain the cross sections at this scale, we must add the terms proportional to ln(m 2 //z 2 ) as 
described in Sec. 1111 A"l in the massless limit and the corresponding terms for the FFN scheme 
contained in Ref. p|. In Fig.^3, we first show the results for the cross section ratios including 
the rescaling due to renormalization and initial-state factorization, i.e. adding Eq. (}2Tj) using 
the coefficients of Eqs. P% |) - (p?Tj) for the massless calculation and the corresponding mass- 
dependent terms of Ref. [2| for the massive calculation. As in Fig. [T^., we show four curves 
corresponding to the massive and massless calculations, in either case without and with 
finite terms subtracted. In this case, the cross section ratios exhibit a much weaker pt 
dependence. The QCD correction (K) factor is somewhat smaller now, but it is still large, 
showing that the perturbative expansion for the gg channel is not converging very well. 
We observe that the massive cross sections converge to the corresponding massless cross 
sections with increasing value of px as in Fig. Ct. The effect of the subtraction of the finite 
terms is slightly smaller, since the added ln(/i 2 /m 2 ) terms apparently have smaller m 2 jp\ 
corrections. The curves for the massive theory lie always above the massless approximation, 
as in Fig. [IJi. 

If it were not for the choice rif = 4 in a s and the fact that the CTEQ6M proton PDFs are 
evolved according to the ZM-VFN scheme, so that the charm and bottom PDFs participate 
in the DGLAP evolution for sufficiently high values of /if, the unsubtracted result a m in 
Fig. \Ih would correspond to the cross section in the genuine FFN scheme. We have seen that 
this theory is characterized by large NLO corrections. This has its origin in the fact that 
the contributions of the would-be collinear divergences related to incoming and outgoing 
charm-quark lines are not yet subtracted, i.e. they are still left at the factorization scales 
y,p — jjilp — m. It is clear that these contributions must be also evaluated at the factorization 
scales fiF = h'f = ^ m T, since our FF for c — > D* + is evolved to this scale and we want to 
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include the contributions from the charm content of the incoming (anti-)proton. To do 
this, we must include the additional contributions proportional to ln(/x|,/m 2 ), ln(fi F /m 2 ) 
and \n(fip/m 2 ) in Eqs. (|35 |1 — (|4(J |1 . The result is shown in Fig. QJ;. The subtracted massive 
calculation (cr^) is represented by the dash-dotted curve and the result of the subtracted 
massless calculation (cr^) by the dotted one. The latter curve again agrees with the ZM- 
VFN result calculated on the basis of Ref. j^J. For comparison, we also show the results 
for o"o and a m from Fig. before the additional logarithmic terms from Eqs. fl35p - (j4*U|l are 
subtracted, and also without subtraction of the finite Acj terms. The lower two curves are 
our final results for the gg channel. The cross section ratios are negative for < 30 GeV 
and rise up to approximately 0.6 at px — 100 GeV. The massive cross section approaches 
its massless limit with increasing value of p?- As we shall see later, the gg channel is only 
important at small values of px\ at higher values of pt, the total cross section is dominated 
by the contribution due to the charm content of the (anti-)proton. Since this dominating 
contribution does not contain finite-m corrections, the effect of mass-dependent terms in the 
gg channel will be suppressed in the final result. 

In Fig. we present results for the qq channel. The cross section ratios shown here 
are normalized to the corresponding LO cross section for the qq initial states. Since this 
normalization differs from the LO gg cross section, the ratios in Figs. ^ and |2 should not 
be added. In Fig. |2t, the ratios with the renormalization and factorization scales equal to 
m are shown. The massless and massive cross sections, i.e. the full (<7o) and the dashed 
(a m ) curves with no subtraction of finite terms and the dotted (cr^) and dashed-dotted (cr^) 
curves with Aq terms subtracted, almost coincide. This means that, for the qq channel, 
the m 2 jp T terms are negligibly small in the considered pt range. We notice that the NLO 
cross section for this channel is negative, except for small values of px- The subtraction 
terms are non-negligible over the whole pt range. In Fig. Eb, the cross section ratio at the 
scale fin = = ^'f = 2wt including the rescaling due to renormalization and initial-state 
factorization according to Eqs. (j2Zj), (JH5|) - (JH%)) is shown. At this scale, the cross section is 
positive for all px values. The influence of the m 2 jp\ terms is somewhat larger now due 
to additional mass- dependent contributions in the terms proportional to ln(/i 2 /m 2 ). Also 
the difference due to the subtraction of the finite Acj terms is larger, since the PDFs and 
the FF are evaluated at much larger scales. We note that the massive ratio is smaller than 
the massless one. If the logarithmic rescaling terms due to internal charm-quark lines in 
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FIG. 2: p(q) + p(q) contribution to pp — > D*^ + X including Abelian and non-Abelian parts, 
normalized to the LO cross section with m = 0. (a): Renormalization and factorization scales are 
fJ-R = [J-F = fJ-'p = Tn (but fixed at 2.1m in a s , PDFs and FF). (b): Same as in part (a), but for 
fJ-R = fJ'F = = ^ m T including rescaling due to renormalization and initial-state factorization, 
(c): Same as in part (a), but including full rescaling to pr = [If = Hp = 2mr. 

Eqs. (jHIj), (JHnj)-((n2J) are subtracted, we obtain the cross section ratios presented in Fig. Efc. 
These results are needed in our final analysis. The dashed-dotted and dotted (oq) 
curves lie very near to each other showing that the m 2 jy% terms are much smaller in this 
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channel. Also in this case, we checked that our result for the massless calculation with 
subtraction and rescaling included agrees perfectly with the ZM-VFN result of Ref. Q- 
For comparison, the unsubtracted massive (cr m ) and massless (<To) results from Fig. |2b are 
again shown in Fig. |2fc. In contrast to the gg channel, the last rescaling due to singular 
internal charm-quark lines leads to an increased ratio. 

Finally, we discuss the m 2 jp\ contributions for the gq and gq channels. For defmiteness, 
we consider the p(q) + p(g) contributions to p + p — > D*^ + X separately for the cases 
where the D*^ meson originates from a charm or an anti-charm quark. We normalize the 
cross sections to the LO gg-channel cross section. Using as scales m in the same way as 
in Fig. Hk, we obtain the full and dashed curves in Fig. corresponding to massless (er ) 
and massive (cr m ) cross sections, respectively. The upper curves are for observed anti-charm 
and the lower ones for observed charm. These two sets of curves differ only little. The 
m 2 power corrections are small over the whole px range. If we rescale to the scale 2m j- 
using Eqs. (|2*7j) . (|H7jl . we obtain the massless (full curves) and massive (dashed curves) 
results for observed charm in Fig. and those for observed anti-charm in Fig. Efc. Here, 
the ratio of cross sections is positive. We observe that the massive cross section is slightly 
larger than the massless one. If we include the full rescaling of logarithmic terms due to 
charm or anti-charm exchange using Eqs. ()27|) and (JMj) with Eqs. (|66| ) - ()7UJI . we obtain the 
lower curves in Fig. |3b and Efc, respectively, for the massive (dashed-dotted curves) and 
massless (dotted curves) cross sections. The latter ones are again in perfect agreement with 
the ZM-VFN results of Ref. jj. These results will be used for the final cross section. 

So far, we have only considered the channels with gluon or light (anti-) quarks in the 
initial state and (anti-)charm quarks in the final state. To these contributions, we must add 
the cross sections with (anti-)charm quarks in the initial state. Such results will be shown 
in the next three figures. From Figs. Eb and c, we see that the gg and gq channels 
yield negative contributions (for px < 30 GeV in the gg case and for all values of pt in the 
gq case), while only the qq channel gives a positive contribution everywhere. Their sum is 
still negative for pt < 65 GeV, as can be seen in Fig. HI where we have plotted the sum 
of these three channels as the solid curve labeled tot. Only when the contributions from 
incoming (anti-)charm quarks are added, we obtain, as a meaningful result, a positive cross 
section, which is shown in Fig. 0] as the dashed-dotted curve labeled tot+ cc. These results 
are for the subtracted massless approach (<7q ) and represent the cross section for inclusive 
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FIG. 3: (a): p(q) + p(g) contribution to pp — > D*^ + X including Abelian and non-Abelian parts, 
normalized to the LO cross section of the gg channel with m = 0. Renormalization and factorization 
scales are hr = [i F = y! F = m (but fixed at 2.1m in a s , PDFs and FF). Lower curves: observed 
charm, upper curves: observed anti-charm, (b): Same as in part (a) for observed charm, but with 
renormalization and factorization scales chosen as (Ir = \ip = fi' F = 2rriT- (c): Same as in part 
(b), but for observed anti-charm. 

jj*± production, including only the contribution due to the fragmentation of a (anti-)charm 
quark, but not the ones due the to fragmentation of a gluon or a light (anti-)quark. In 
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FIG. 4: Partial cross sections for pp — ► D + X m. the ZM-VFN scheme. Renormalization and 
factorization scales are = [af = = 2my. The sum of the gg, qq, gq and gq channels is 
labeled tot, the total cross section including incoming (anti-)charm quarks tot + cc. 

addition to these results, we also show in Fig. 0] the partial cross sections for the gg, qq and 
g(q + q) channels. The latter channel contains all contributions with a gluon and a light 
(anti-)quark coming from the (anti-)proton. We see that the sum of all components yields 
a cross section with a smooth, steeply falling px dependence. The comparison with the 
equivalent results with massive (anti-)charm quarks is made in Fig. where we show for 
each partial result the ratio with respect to the full result including the contribution from 
cc initial states. The upper and lower curves correspond to the subtracted massive (cr^) 



28 



_Q 4 l~. . I .... I .... I . I .. I .... I .... I . I . I I ... I I I ... I .... 1 

10 20 30 40 50 60 70 80 90 100 

p T [GeV] 

FIG. 5: Partial cross sections normalized to the result for all contributions including the cc initial 
state. Lower and upper curves correspond to the massless and massive calculations, respectively. 

and massless (cq) calculations, respectively. As we can see, the w? jp\ corrections in the 
sum of all contributions with gluons or light (anti-) quarks in the initial state labeled tot are 
modest, except at small values of px- At the smallest pt value considered, pt = 5 GeV, the 
massless cross section is only increased by roughly 20% when these corrections are included, 
although the contribution labeled tot is substantially increased, from —773 pb GeV -1 for the 
massless calculation to —386 pb GeV -1 for the massive one. We emphasize that all results 
in Figs. |U and El refer to the scale choice /ir = Hf = ^'f = ^ m T- 
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IV. COMPARISON WITH CDF DATA 



Before we can compare our final results with experimental data, we have to add another 
contribution, which was not yet discussed, but is non-negligible for the experimental situa- 
tion at the Tevatron: in fact, the observed D*^ meson may also be produced through the 
fragmentation of a gluon or a light (anti-)quark. Appropriate g,q,q — > D*^ FFs are con- 
tained in the OPAL set of Ref. |3|, where they are generated via NLO evolution assuming 
that they vanish at the starting scale. In the ZM-VFN scheme, gluon and light (anti-)quark 
fragmentation already contributes at LO. By contrast, these types of contribution do not 
exist in the FFN scheme, where the collinear g — > cc splitting is treated perturbatively 
instead, starting only in NLO. Since the GM-VFN scheme is to be constructed in such a 
way that it merges with the ZM-VFN scheme at large values of pt, it is clear that the gluon 
and light (anti-) quark fragmentation contributions must be accommodated in the GM-VFN 
framework as well. A certain part of these contributions carries mass dependence, because 
of internal charm-quark lines and external ones that do not lead to fragmentation. In the 
FFN scheme, the analogous mass dependence only comes in beyond NLO. Therefore, we 
ignore this mass dependence for the time being and adopt the gluon and light (anti-)quark 
fragmentation contributions from the ZM-VFN analysis [20]. These contributions amount 
to slightly more than 30%, almost independent of pt, as can be seen in Fig. |U1 The bulk is 
due to gluon fragmentation. For photoproduction in ep and 77 collisions, this contribution 



was found to be negligible (la, LLZI, • 



The effect of mass- dependent terms is very much reduced in the final cross section, since 
the parts which have to be calculated with m = are large. Therefore, one cannot expect 
that mass terms would increase the theoretical predictions towards cross sections as high as 
observed in the CDF experiment [l| . The size of the mass- dependent terms is visualized 
in Fig. [?l where we show the ratios of cross sections calculated with m 7^ to those with 
m — 0. For the (negative) contributions due to incoming gluons and light (anti-) quarks, 
mass- dependent terms lead to a reduction in size by 50% at pt = 5 GeV. At this value of pt, 
the ratio of massive over massless results is reduced to 1.19 and 1.13 when the contributions 
from incoming (anti-)charm quarks and from g — > D* ± fragmentation, respectively, are 
included. 

Since the effect of mass-dependent terms is very much reduced in the final cross sec- 
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FIG. 6: px spectrum for pp — > D* ± + X including all contributions. Renormalization and factor- 
ization scales are [Xr = \ip = fx' F = 2mr- Dash-dotted lines: the (negative!) contributions with 
gluons and light (anti-)quarks in the initial state; dashed lines: including (anti-)charm quarks in 
the initial state; full lines: including the contribution from g,q,q — > D*^ fragmentation. Upper 
and lower curves correspond to the massive and massless calculations, respectively. 

tion, it is clear that a variation of the value of the charm-quark mass does not contribute 
much to the theoretical uncertainty. Whereas the sum of the contributions due to gluons 
and light (anti-)quarks in the initial state varies by —20.5% (+17.8%) at p T = 5 GeV if 
m is changed from our default value of 1.5 GeV to 1.2 (1.8) GeV, the cross section also 
including those terms which can be calculated in the massless approach only (i.e. those with 
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FIG. 7: Ratios of the subtracted massive and massless differential cross sections da/px for 
pp — ► Z)*^ + X with |y| < 1 including (a) all contributions (full line), (b) all contributions with in- 
coming gluons, light (anti-)quarks and (anti-)charm quarks, but without the g,q,q —* D*^ fragmen- 
tation contributions (dashed line), and (c) the sum of the contributions with only gluons and light 
(anti-)quarks in the initial state (dash-dotted line). Renormalization and factorization scales are 
Hr = = Hp = 2m T - 

(anti-)charm quarks in the initial state and from g,q,q — > -D* ± fragmentation) varies by 
only —2.5% (+1.8%). The theoretical uncertainty is, therefore, dominated by the choice of 
the renormalization and factorization scales. The results presented in Fig. |H1 were obtained 
for = fx p = fi' F = 2rriT- A conservative mode of scale variation frequently encountered 
in the literature is to independently vary the values of hr/ttiti lif/wit and n' F / m T between 
1/2 and 2 about the default value 1. We will adopt this prescription for the comparison 
with experimental data from the CDF Collaboration This leads to large changes at 
Pt = 5 GeV (for the massless calculation in the range between —41% and +30%, and for 
the massive calculation between —46% and +56%). At px — 25 GeV, the variations are 
smaller and cover the range ±20%. Minimal values for da/dpT are obtained for maximal 
values of lir and li' f and minimal value of tip, and maximum values are reached for minimal 
values of hr and fjf F and maximal value of Hp. These results are presented in Fig. |H] as a 
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FIG. 8: Variation of the px spectrum of pp — > D* ± + X with the renormalization and factorization 
scales in the GM-VFN scheme. The central solid curve is for /j>r = hf = n' F = mj-, the upper 
and lower dashed curves represent the maximum and minimum cross sections found by varying 
fj>R, [ip and [i' F independently within a factor of 2 up and down relative to the central values. The 
contributions from (anti-)charm quarks in the initial state and from g,q,q — * -D* ± fragmentation 
are included. The data points are from CDF |lj for the average of D* + and D*~ production. 

band of predicted values. A comparison with experimental data from CDF also shown in 
this figure, shows that the data prefer a small renormalization scale and a large initial-state 
factorization scale. However, even for the choice fiR = rriT/2 and fi F = 2my, the theoretical 
results in the GM-VFN scheme still undershoot the measured cross sections by 20% to 34%. 
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Residual sources of theoretical uncertainty include the variations in the available PDF 
and FF sets. The new generation of NLO proton PDF sets exhibit only minor differences. 
Similarly, the two NLO FF sets that were determined in Ref. [^| by fitting slightly incompat- 
ible ALEPH and OPAL data yield rather similar predictions, the difference being of order 
10% or less in the kinematic range of the CDF data The inclusion of these additional 
errors would broaden the theoretical error band in Fig. |H] only insignificantly. 

We stress that our results do not include a contribution due to bottom-quark production. 
This part can be identified in the experimental analysis and was, in fact, removed from the 
CDF data shown in Fig. |H1 



V. SUMMARY AND CONCLUSIONS 



In this work, we calculated the NLO corrections to inclusive charm production in pp col- 
lisions in two approaches using massive or massless charm quarks. By deriving the massless 
limit from the massive theory (FFN scheme), we could show that this limit differs from 
the genuine massless version with MS factorization (ZM-VFN scheme) by finite corrections. 
We adjusted subtraction terms and thus managed to establish a massive theory with MS 
subtraction (GM-VFN scheme) which approaches the massless theory (ZM-VFN scheme) 
with increasing transverse momentum. With these results and including the contributions 
where a (anti-)charm quark occurs as an incoming parton and those where a gluon or light 
(anti-)quark fragments, we calculated the inclusive D* ± cross section in pp collisions at 
= 1.96 TeV using realistic non-perturbative FFs, which are subject to proper DGLAP 
evolution and manifestly universal ^ . Our central prediction somewhat undershoots a 
recent measurement by the CDF Collaboration [lj , but reasonable agreement can be reached 
by adjusting the renormalization and factorization scales within a plausible range of toler- 
ance. 

We made the observation that, in contrast to other experimental situations such as 77 or 
7P collisions, the contribution from g — > D*^ fragmentation is large. The FFs which we used 
here had been obtained from fits to data from the CERN Large Electron-Positron Collider 
(LEP). One may speculate, therefore, whether these data leave enough room to adjust the 
g — > D*^ FF and improve the agreement with the Tevatron data. We plan to come back to 
this question in a future publication. 
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APPENDIX A: COEFFICIENTS FOR SUBPROCESS g + g -»■ c + X 



In this appendix, we list the coefficients c, needed for the calculation of the cross section 
for the inclusive production of charm in gg co 
m — > of the cross sections calculated in Ref. 



lisions. They are obtained by taking the limit 
~\. These limits are compared with the results 
from Ref. [20j in order to obtain the subtraction terms Aq. To shorten the expressions, we 
make use of the abbreviations 

X = 1 — vw , Y = 1 — v + vw , Vi = i — v , (Al) 

and write down only non-zero contributions. The coefficients given below have to be used 
together with Eqs. (jl2j) and ([20)): their colour decomposition is defined in Eq. (j!6|) . We start 
with the Abelian coefficients c? ed . 
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We note that these results for the Abelian part agree with the coefficients for 7+7 — > c+X 
in Ref. [lfij . if the normalization factor C(s) in Ref. Q| is replaced by C(s) = 1. In that 
work, the zero-mass limit was derived from the FFN cross sections of Ref. 
with the ZM-VFN ones of Ref. bfll. 
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2. c° q Coefficients 
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4v (1 — 5v)w 



v\ 



2v 2 {7 + v) 8(1 - 2v + 2v 2 ) 2(2 - 3v + 2v 2 + v 3 + 2v 4 ) 
v\ V\VW 2 v\vw 

2(-9 + v)v 2 w 4v(2-2v + v 2 ) 
+ v 2 + v^X ' 

-2(2 - 8v + 4v 2 - 9v 3 + 3v 4 ) 8(1 - 2v + 2v 2 ) 



V 2 V V\VW 2 

2(-4 + 2v + 3v 2 - 8v 3 + 3v 4 ) 2(2 - 8v + 8v 2 - llv 3 + v 4 )w 

^ 2 ^ 2 

vfvw vfv 
4v x v 4v{-3 + 2v) 4v{-3 + v 2 ) 
+ + x2 

4u?u 4wiw 2 2v(3-6v + 4t; 2 ) kq 

-10 + 31u - 22v 2 + 7w 3 + 2v 4 _ 8 2(1 - Qv + Qv 2 ) 

vfv w 2 vw 

2(2-5v + v 2 )(3-2v + v 2 )w \2v x v 2v{-\2 + 7v) v{-3 - 9v + 10v 2 



vfv X 3 X 2 Vl X 

8v\v 4uiu(-l + 3u) _ u(-19 + 6v 2 ) kq 

~ ~y3~ ^ y 2 y ^ 11 ' 

2 10 ' 

2(2 - 8v + 10w 2 - 13v 3 + w 4 ) 8(1 - 2v + 2v 2 ) _ 2(-5 + \2v - 13w 2 + 2v 3 ) 
vfv v\vw 2 v\w 

4(-l + 4v - 4v 2 + 5v 3 )w 4v x v 4v{-3 + 2v) 4v(-3 + v 2 ) 
+ v\v ~X~ 3 x 2 + v^X 

4v\v 4v x v 2 2v{3 -Qv + 4v 2 ) 



y3 y 2 Y 
8_ 8 

Vi v 



6 — lQv + lQv 2 , 



4 4 
10 - 2w + , 

V\ V 

4 4 
8 + 2v + - . 

V\ V 
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4. Coefficient for the Quark-Loop Contribution 



cf = Acf , 

Acf = --v(l - v) . (A5) 

APPENDIX B: COEFFICIENTS FOR SUBPROCESS q + q^ c + c 

In the following, we list the coefficients needed for the calculation of the cross section for 
the inclusive production of charm in qq collisions. As before, the coefficients are obtained 
from the cross sections of Ref. I7f by taking the limit m — > 0, whereas the subtraction terms 

n 

Aq are deduced by a comparison with the results of Ref. [2fJ. The coefficients given in this 
appendix determine the cross section according to Eqs. (fT2|) and PU|) and using the color 
decomposition of Eq. (|49J1 . We start with the Abelian coefficients cf. 

1. c| f Coefficients 



cf = -[2Av\nv 1 + 6(3-Qv + 2v 2 )\n 2 v 1 + 5(-15 + 2n 2 )(l-2v + 2v 2 ) 

- 24(3 - 6v + 7v 2 ) \n 2 v + 3(— 11 + Uv - 6v 2 + 28(1 - 2v + 2v 2 ) \nv-y) Inv] 
+ Acf , 

Acf = +\nv + \n 2 v)T q (v) , 

cf = (-9 - 8lnv + 4lnvx)r q (v) , 

cf = (-3-20lnv + 24:\nv 1 )T q (v) + Acf , 
Acf = -4(l + 21nv)r ? (v), 

£f = -12t», 

cf = 20r fl (t;) + Acf , 
Acf = -8t» , 

cf _ 2(-l + 34w-39w 2 + 8t; 3 ) 4(1 - 2w + 2v 2 ) 8v 2 w 8v 3 w 2 2v 
Vi w Vi Vi X 

Avjv Aviv 2 2(55w-74v 2 + 28v 3 ) f 
+ ~Y~z y2 y ^ 5 ' 
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8vv 2 8v 2 v x A(3v - 6v 2 + Av 3 ) 

Av -I — - 

Y 3 Y 2 Y ' 

-Av(7 - lOu + 4w 2 ) 8(2 - v)v 2 w 8v 3 w 2 



+ 

Vl V! Vi 

8vvf 8v 2 vi A(v - Qv 2 + Av 3 ) 

- 20v + ~ 3 yT + y ' 

Av(-3 + 2v + 2v 2 ) 8v 2 w 8v 3 w 2 16(5u - 8v 2 + Av 3 ) 

v\ v\ vi Y 

4(1 - 2v + 2v 2 ) 9 16(5v -8v 2 + Av 3 ) 
AOv + — + 8v 2 w - v 



w Y 
2(-l + 22v - 31v 2 + 12v 3 ) 4(1 - 2v + 2v 2 ) 8(2 - v)v 2 w 8v 3 w 2 2v 

vi w vi vi X 

Avv\ 2(-15u + 10w 2 + 4w 3 ) f 

H 775 77^ I 77 I" Ac ln i 



Ar cf 

-10 + 9v + 3v 2 2 6v 8vv 2 12(3v - Av 2 + v 3 ) 2(-llv + 3v 3 ) 



+ Ac cf 



vi w X Y 3 Y 2 Y 

li ) 



-Ar ci 
2 10 ' 

2(1 - lOu + 15w 2 - 8v 3 ) 4(1 - 2w + 2w 2 ) 8v 2 w 8v 3 w 2 2v 
vi w vi vi X 

Av\v Aviv 2 2(v - 6v 2 + Av 3 ) 

8(-l + 2w + 2w 2 ), 
-8(3 - 6v + 8v 2 ) , 

lQvj. (Bl) 



Coefficients 



- [ - 18u In vi + 9(-l + 2v) In 2 Vl - 2(-85 + 9tt 2 ) (1 - 2v + 2v 2 ) 
9 

+ 18(7- Uv + 16v 2 ) In 2 u - 36(-l +v+ (3 - 6u + 6u 2 ) lnui) In v] , 
22 r\ 

16r g (v) Inv , 

0, 

0, 
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10 



ca 
c ll 



~ca 
c ll 



12 



13 



11 



-Av (4 + f ) - 4w 2 u; - 



4wf 4(3f - 5v 2 + 2w 3 ) 2(19u - 28v 2 + 12v 3 



Y 3 



Y 2 



Y 



Av , 

8w? 8(3w - 5v 2 + 2v 3 ) A(9v - 12v 2 + Av' 3 
I8v - + -i— 



r 3 

-2v(-7 + 2v) -4v 2 u; 
-4t> (4 + v) - Av 2 w + 



Y 2 Y 
2(17v - 26v 2 + 12v 3 ) 



Y 

2(17v - 26w 2 + 12v 3 
Y 



Avv\ A(3v - 5v 2 + 2v 3 ) Avv x 

—OV TT^T H 777^ h 



Y 3 



Y 2 



2(2 - 5v + 2v 2 ) + Av 2 w - 



Y ' 

2w 12vvf \2vv\ A(Av - 7v 2 + 2v 



X 



Y 3 



Y 2 



+ 



Y 



9 Aw 2 A(3v - 5v 2 + 2v 3 ) 2(9v - 12v 2 + Av 2 
4 -2 + v)v + Av 2 w + — - - — — : '- + 



Y 3 



Y 2 



Y 



-2 + Av- 16v' 
4(3 - 6v + 8v 2 ) 
-2 + Av. 



(B2) 



3. Coefficient for the Quark-Loop Contribution 



20 



-n f r q {v) 



A 



(B3) 



APPENDIX C: COEFFICIENTS FOR SUBPROCESS g + q -»■ c + c + q 

This appendix contains the coefficients needed for the calculation of the cross section for 
g + q — > c + c + g with an observed charm quark in the final state according to Eq. (|12j). 
We note that there are no subtraction terms for the gq channel. In the following, again 
only non-zero coefficients are written down. According to the color decomposition defined 
in Eq. (fo3|). we present the Abelian and non-Abelian coefficients cf and cf 1 separately. 
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1. c| Coefficients 



c 



-2(1 - 2v + Qv 2 - 5v 3 + 2v 4 ) 2(1 - 2v + Av 2 - 3v 3 + v 4 ) 
5 vf vfw 

2(1 - 2v + 6v 2 - 5v 3 + 2v 4 )w v -3v + v 2 v x v , 3v - 2v 2 

I O I ~ Trn I 



cf 



10 



v 2 2X 2 2v x X 2Y 2 2Y 



2v(l + v 2 ) -1-v 2 Av 3 w 



r rf _ ^_ L _|_ 

cf — 2t> 4t> 2 u> ««i 3t> — 2v 2 

c 7 = ' TP! ^ t7 ' 

V\ Vi Y Y 

cf -2v(3 - 3v + 2v 2 ) 1 + v 2 Av(2-3v + 2v 2 )w 

C 8 2 ' 2 ' 2 



-4(1 - 2v + 5v 2 - Av 3 + v 4 ) 2(1 - 2v + Av 2 - 3v 3 + v 4 



9 — 9 ' 9 

v{ v{w 
Av 2 (5 — 5v + v 2 )w 16v 2 w 2 



+ 



-2(-l + 2f - v 2 + 2v 4 ) 2(1 -2v + Av 2 - 3v 3 + v 4 ) 



vi viw 

2(1 - 2v - Av 2 + 5f 3 + 2v 4 )w lQv 2 w 2 v -3v + v 2 
-I — * '- 1 1 1 

vj v x 2X 2 2v x X 

V\V 3v — 2v 2 
+ 



2F 2 2Y ' 
2v 2 (A-5v + 2v 2 ) l + 2v-2v 2 



f 2 2-u; 
-2 + Av + 7v 2 - llv 3 + Av 4 )w v 5v - 4v 2 t^u ( -8u + 3v 2 



v 2 2X 2 2^X F 2 2F 

„ cf 2(1 - 2f + 6v 2 - 5f 3 + 2f 4 ) 2(1 - 2v + Av 2 - 3v 3 + v 4 ) 

C ll ~~ 2 2 

V i viw 
2(1 -2v + Qv 2 - 5v 3 + 2v 4 )w v 3v - v 2 v x v -3v + 2v 2 

v 2 2X 2 + 2 Vl X + 2~Y 2+ 2~Y ' ( ' 

We have checked that the following relations between these results for the Cp part and the 
coefficients for the subprocess 7 + g — > c + c + q in Ref. ^1 are satisfied, if the normalization 
factors C(s), C q (s), C cq (s) and Cp are replaced by unity: cf = cf 1 +cf 2 —2cf 3 . Furthermore, 
our results in Appendix ID II fulfill the relations cf ,c ~ c = Acf 3 and cf + cf' c ~ c /2 = cf 1 + cf 2 . 

The decomposition of the form e 2 Qi + e 2 Q2 — e c e q Q3 for jq — > Q turns into one of the 
form Qi + Q 2 — 2Q 3 for gq — > Q. The interference part Q3 is anti-symmetric and changes 
sign in the analogous subprocesses with an anti-quark q, i.e. we have e 2 Qi + e 2 Q2 + Zc^qQz 



AA 



for 7g — > Q and Qi + Q 2 + 2Q S for gg — > Q. 



c*f a Coefficients 



v 2 (2 — 2w + w 2 ) 4t> 2 (l — + -y 2 )u> 6f 4 w 2 4w 4 w 3 u u 

5 t; 2 + i; 2 v\ v\ 2X + 2Y" 

ca f(l + f) v 2 (5 — v + 2v 2 )w Av 3 (l + v)w 2 Av A w 3 

vf vf vf vf 

V 2 W V 

cf = -v + - , 



Vl Y 



Co r, + r> r, + 



v (-1 + 4f - 3v 2 + v 3 ) 2v (-1 + 4v - 3v 2 + 2v 3 )w 6v 4 w 2 4t> 4 w 3 



r ca — ; __ _i_ 



2 1 2 2 1 2 ' 

i>f ff vf vf 

-l + 2v-v 2 + v 4 v 2 (-l+v + 2v 2 )w 2v 2 (-2 + 2v + v 2 )w 2 



V 2 V 2 



u l u l "1 

l-v + 5v 2 v 2 (ll - lit) + 2v 2 )w Av 2 (l-v + v 2 )w 2 4v 4 w 3 

u 10 "T" 2 2 2 

V\ vf vf vf 



V V 

2X + 2Y 



, 2(-l + v + t> 3 W 2v 3 (-l + 2vW 2 2w 4 w 3 w w 

C 1 = 1 + ^ + — + g " — -- + —, 

vi vf vf X 2Y 

~ ca v 2 (2 — 2v + v 2 ) Av 2 (l — v + v 2 )w Qv 4 w 2 4v 4 w 3 v v 
Cn= v\ v\ ^r + 2X _ 2Y- (C2) 

APPENDIX D: COEFFICIENTS FOR SUBPROCESS g + q -> c + c + q 

The cross sections for the subprocesses g + g — > c + c + g and g + q — > c + c + q, where the 
quark in the initial state is replaced by an anti-quark, are related, but not identical. The 
differences between the corresponding coefficients, c^~ c , will be presented below. They have 
to be combined with the coefficients for g + g — > c + c + q given in Appendix to give those 
for g + q — > c + c + q, according to c^{g + q — > c + c + q) = Ci(g + q — > c + c + q) + c^ _c , to be 
inserted in Eqs. f!12|) and (|65j) . Again, there are no subtraction terms, and we present only 
the non-zero coefficients. Note, in particular, that for both colour factors cfZ^ and C12-14 
are zero. 
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1. cf Coefficients 



2. cf" Coefficients 



4v 2 8v 2 w 



cf, c—c 

C 6 = 

Vl V! 

cf, c—c cf, c—c 

c 7 — C 6 ' 



Cg f ' c c = 8f — 16fU7 , 

cfjC -_ c _ 4(2 - 2t; + 3t; 2 ) 24t; 2 w 32v 2 w 2 

c 9 ' ~~ I ) 

Vl Vi Vi 

cfi5 _ c -4(2 - 2v + 5v 2 ) 40v 2 w 32v 2 w 2 
c io ~~ I ' 

rf r-r 4t> 4t> _ , 

cf{ c = 8 - 16™ - — + — . (Dl) 



-3v 2 3v 2 w 
+ 



2vi 



Vl 



= c R 



= —3v + 6vw , 



'10 



ca,c— c 
C ll 



3(2 - 2w + 3t> 2 ) 9u 2 w 



2«i 

3(2 - 2v + 5v 2 



12v 2 u> 2 



2ui 
-3 + 6w + 



15v 2 w 12w 2 w 2 
+ 



Vl 



Vl 



3v 
2X 



3v 
2Y 



(D2) 
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